arXiv:1505.05036vl [astro-ph.IM] 19 May 2015 


A four-pole power-combiner design for far-infrared and 
submillimeter spectroscopy 

Giuseppe Cataldo a,b ’*, Samuel H. Moseley a , Edward J. Wollack a 

a NASA Goddard Space Flight Center, United States of America 
b Massachusetts Institute of Technology, United States of America 


Abstract 

The far-infrared and submillimeter portions of the electromagnetic spectrum 
provide a unique view of the astrophysical processes present in the early uni¬ 
verse. Micro-Spec (p-Spec), a high-efficiency direct-detection spectrometer con¬ 
cept working in the 450-1000-pm wavelength range, will enable a wide range 
of spaceflight missions that would otherwise be challenging due to the large 
size of current instruments and the required spectral resolution and sensitivity. 
This paper focuses on the ft -Spec two-dimensional multimode region, where the 
light of different wavelengths diffracts and converges onto a set of detectors. 
A two-step optimization process is used to generate geometrical configurations 
given specific requirements on spectrometer size, operating spectral range, and 
performance. The canonically employed focal-plane constraints for the power 
combiner were removed to probe the design space in its entirety. A new four- 
stigmatic-point optical design solution is identified and explored for use in far- 
infrared and submillimeter spectroscopy. 


1. Introduction 

Far-infrared (IR) and submillimeter (15 pm to 1 mm) spectroscopy pro¬ 
vides a powerful tool to probe a wide range of environments in the universe. 
In the past thirty years, discoveries made by several space-based observatories 
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have provided unique insights into physical processes leading to the evolution 
of the universe and its contents. This information is encoded in a variety of 
molecular and fine structure lines; observations of such spectral lines enable 
the exploration of galaxies at high redshifts. The fine structure lines of abun¬ 
dant elements (carbon, nitrogen, and oxygen), for example, allow tracing the 
obscured star formation and Active Galactic Nuclei (AGN) activity into the 
high-redshift universe. One can measure galaxy redshifts and determine their 
elemental abundances and physical conditions out to redshifts z > 5. In spite of 
this, a number of questions remain unanswered regarding the very early steps 
of the universe as well as galactic, stellar, and planetary formation. The ability 
to explore this rich spectral region has been limited by the size and cost of the 
cryogenic spectrometers required to carry out these measurements. The work 
proposed here specifically addresses the need for integrated spectrometers and 
background-limited far-IR direct detectors. For space-borne astrophysics sys¬ 
tems, the specific requirements are shown in Table |T| and compared against the 
current state of the art [l]. 


Table 1: Summary of far-IR cryogenic spectrometer and detector array require¬ 
ments and comparison with current state of the art p]. 


Metric 

State of the art 

Requirements 

Wavelength, A 

250 - 700 pm 

220 - 2000 pm 

Noise Equivalent Power, NEP 

10" 19 W /VBz 

< 10" 20 W/VHz 

Spectral resolution, 1Z 

> 100 

> 1200 

Detective Quantum Efficiency, DQE 

~ 15% 

> 90% 

Time constant, r 

100 ms 

< 10 ms 


In order to realize the goals outlined in Table [T] a high-performance inte¬ 
grated spectrometer module, Micro-Spec (p-Spec), operating in the 450-1000- 
pm (300-650-GHz) range is proposed. p-Spec can be compared to a grating 
spectrometer [21 Sj .41151 IS] 7 in which a plane wave is reflected from the grating and 
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the phase of each partial wave scattered from the rulings is a linear function of 
position across the grating. An example of planar Rowland grating architecture 
is Z-Spec, in which propagation occurs in parallel-plate waveguides [HIEl SEED]. 
Another comparison can be made with one-dimensional bootlace lenses found in 
microwave practice pum ns nu, which /Lt-Spec builds upon for submillimeter 
wave applications. Finally, another variation used at millimeter wavelengths, 
which does not rely on optical interference as in grating spectrometers, is a 
narrow-band filter-bank spectrometer. Examples realized in superconducting 
transmission lines are SuperSpec H5JHIE1IZ] and the Delft SRON High-Z Map¬ 
per (DESHIMA) [IS]. 

/z-Spec differs from these approaches by the order of processing of the light 
in the spectrometer. In /i-Spec (Fig. [I]), the incoming radiation collected by the 
telescope is coupled to the spectrometer via a broadband dual-slot antenna used 
in conjunction with a hyperhemispherical silicon lens and directed to a series 
of power splitters and a delay network made of superconducting microstrip 
transmission lines. Analogous to the Rowland grating [2j, the delay network 
creates a phase retardation across the input to a planar-waveguide multimode 
region, which has two internal planar antenna arrays, one for transmitting and 
one for receiving the radiation as a function of wavelength. Absorber structures 
lining the multimode region terminate the power emitted into large angles or 
reflected from the receiver antenna array. An array of planar feed structures is 
employed to couple the radiation to the multimode region and concentrates the 
power along the focal surface with different wavelengths at different locations. 
The outputs are connected to a bank of order-sorting filters which terminate 
the power in an array of microwave kinetic inductance detectors (MKIDs) for 
detection and read-out. The entire spectrometer circuit is integrated ona~ 10- 
cm 2 silicon chip (i.e., the hyperhemispherical lens, relay optics, and telescope are 
not on the chip and are part of the instrument system). This compact footprint 
is accomplished through the use of single-mode microstrip delay lines, which 
can compactly be folded on the silicon wafer and reduce the required physical 
line length by a factor of the medium’s effective refractive index. 
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Figure 1: Layout of an individual /x-Spec wafer. The radiation is coupled into 
the instrument through a broadband antenna and is then transmitted through 
a superconducting transmission line to a divider and a phase delay network. 
The spectrum enters the multimode region through an array of feeds which 
concentrates the power along the focal surface with different wavelengths at 
different locations. The receivers are connected to a bank of order-sorting filters 
and MKID detectors m- Multiple spectrometer wafers can be packaged and 
potentially used in defining an instrument system. 
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The frequency range of the implementation presented here is limited to wave¬ 
lengths A > 250 pm by the gap frequency of currently available low-loss super¬ 
conductors. These include niobium (Nb) and niobium-titanium nitride (NbTiN) 
for the transmission line structures, and molybdenum nitride (MoN) for the de¬ 
tectors. This paper will describe the design process of the p-Spec multimode 
region and illustrate the results in terms of geometry, imaging quality, and effi¬ 
ciency. 

2. Multimode region design 

In designing a spectrometer, it is possible to define points on the focal plane 
where the phase error of the diffracted light is identically equal to zero. These 
points are called stigmatic points. Increasing the number of such points on 
the focal surface presents the advantage of improving a spectrometer’s imaging 
quality, which results in lowering the overall phase error on the entire focal plane 
and increasing the usable spectral bandwidth. As a consequence, the number 
of spectrometer channels and the resolving power can grow. 

Examples of designs with two stigmatic points can be found in the litera¬ 
ture BUS El- We built upon these designs to generate spectrometer concepts 
with an increased number of stigmatic points. A three-stigmatic-point prototype 
version with a resolving power TZ = 65 in first order (M = 1) was designed IIS] 
and built, and is currently under evaluation at the NASA Goddard Space Flight 
Center. Additional designs are described in [2U] for configurations with resolv¬ 
ing powers equal to TZ = 260 and TZ = 520 in higher order (M > 1). These 
designs were obtained through a constrained optimization process in which zero 
phase error was imposed on three preselected points. However, a fourth stig¬ 
matic point was observed beyond the angular range in use. It is the purpose of 
the work presented in this paper to show how to use this additional degree of 
freedom to increase the number of spectrometer channels and resolving power. 
In the following section, therefore, we describe a design for TZ ~ 260 in first 
order with four stigmatic points, as a result of an unconstrained optimization 
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process which maximizes the instrument resolving power and minimizes the 
root-mean-square (RMS) phase error on the focal plane. 

2.1. Problem formulation 

As explained in jT5j , the design variables are the x and y coordinates of the 
N e emitters’ centers and the electrical path lengths in silicon, Rf, for each feed’s 
electrical delay (Fig. [2j. The resolving power is defined as 

U = M-N e , (1) 

where M is the order of the spectrometer and N e is a power of 2, given the 
structure of the power divider network (Fig. [I]). 

The first step of the design consists of finding the maximum achievable 
resolving power, lZ m axi as a function of M and the relative emitter pitch, 
r ] = p/Xavg (P = emitter pitch, X avg = central wavelength associated with 
the geometric average frequency), given specific requirements on spectrometer 
radius, R, and operating spectral range, as well as certain constraints on per¬ 
formance. 

The formulation of this mixed integer non-linear problem is as follows: 


max lZmax (M, 7]) 

= 

M/p- R/X avg 

(2) 

subject to H e (M,ij) 

< 

R 

(3) 

H r (M, rj) 

< 

7T R 

(4) 

llmax(M,p) 

> 

200 

(5) 

V 

> 

0 

(6) 

M 

> 

1, M integer 

(7) 


Eq. © imposes that the width of the emitter array, H ei be less than or equal to 
the spectrometer radius so that there is no aberration; Eq. Q lets the receiver 
array be as large as the focal plane to maximize its utilization; and Eq. ([5]) sets 
a minimum value for the required maximum resolving power, lZ m ax , thereby 
eliminating solutions in which we are not interested. Finally we note that, to 
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Figure 2: Simplified representation of the grating geometry. On the left side 
three radiators can be seen, which point to the blaze point, B. The solid lines 
represent the paths to four stigmatic points from a generic z-th radiator’s phase 
center (Ru, i? 2 i, • ■ •) and the central reference feed (i?i, i? 2 , • ■ •)■ 

enlarge the tradespace, the constraint of equal emitter and receiver pitch, used 
in our previous designs [HJ |2Q], was removed. The entire problem was solved 
with a Branch and Bound algorithm [2l| using the Interior Point OPTimizer 
(IPOPT) [22] and Coin-or Branch and Cut (CBC) [22i solvers. 

Table [2] shows the requirements on spectrometer size and spectral range 
used for this problem. The minimum and maximum frequencies are no longer 
associated with any stigmatic point and they fall within the spectral range 
defined in Table [T] The average frequency was computed as their geometric 
mean. 

The objective spaces as a function of M and rj are shown in Fig. [3] On 
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the left (Fig. 3a), it is possible to visualize the feasible objective space of the 
optimization problem described above for lZ ma x along with the active constraint, 
Eq. 0 (blue area). The feasible solutions of Eq. (0 populate that part of 
the contour plot above the blue area, whereas the optimal solutions lie at the 
intersection of the active constraint with the largest TZ max ss 275, one for each 
value of M. One can thus choose the highest desired spectrometer order, M max , 
which is associated with the highest design frequency band. The frequency 
bands corresponding to orders M < M rnax can be calculated by scaling the 
highest one by a factor of M/M max . We decided to investigate the first-order 
case as an example of simple and robust system. Higher-resolution instruments 
will certainly require higher-order operations. Table [3] shows the values of the 
design variables associated with this particular optimal solution as well as the 
values of Eqs. @-0, which satisfy these constraints. 

After solving this problem, it was possible to compute the number of emit¬ 
ters, N e = TZmax/M, and round it down to a power of 2. According to Eq. Q, 
this causes the actual resolving power, 1Z , calculated with this updated value 
of N e . to be lower than lZ max ■ The plot in Fig. 3b shows the values of 7 Z (red 
line) and the values of M and rj that would make such realizations possible. In 
particular, for a first-order design (M = 1 ), the optimal solution corresponds to 
a resolving power 7 Z = 257 with a relative emitter diameter 77 = 0.2916. 


Table 2: Requirements on spectrometer size and spectral range. 


Parameter 

Symbol 

Unit 

Value 

Multimode region radius 

R 

cm 

1.25 

Minimum frequency 

fmin 

GHz 

430.0 

Maximum frequency 

fmax 

GHz 

650.0 

Average frequency 

favg 

GHz 

528.7 


The second step toward determining the optimal solution in terms of the 
above-mentioned design variables consists of minimizing the overall RMS phase 








M 


(a) The contour plot depicts all the fea¬ 
sible values of IZmax associated with dif¬ 
ferent values of M and rj. For M = 1 
and 7/ = 0.2916, the optimal solution is 
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(b) The contour plot represents the fea¬ 
sible values of 7 Z for all orders given as 
powers of 2. The optimal solution for 
M = 1 and 77 = 0.2916 is associated with 
a resolving power 1Z = 257. 

Figure 3: Objective spaces of Problem ([ 2 ]). The blue area represents the infea¬ 
sible region corresponding to the active constraint in Eq. Q. Both plots show 
that, for each order M, several feasible solutions exist for different values of 77. 
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Table 3: Optimal solution of Eq. ([2]) for M = 1. 


Variable 

Symbol 

Unit 

Value 

Spectrometer order 

M 

- 

1 

Relative emitter diameter 

V 

- 

0.2916 

Maximum resolving power 

max 

- 

275 

Emitter array width 

H e 

cm 

1.2452 

Receiver array width 

H r 

cm 

3.5674 

Emitter pitch 

P 

cm 

0.0048 

Receiver pitch 

s 

cm 

0.0167 

Number of emitters 

N e 

- 

257 

Number of receivers 

N r 

- 

149 

Resolving power 

n 

- 

257 


error, </?rmS; through a figure of merit representing the area subtended by </?rms 
on the focal plane as follows: 


min / ^rms(/ 3) d/3, 


with 


T’RMS = 


N e 

E 


[ipi{xi,yi,Rf,/3) - M/3))] 2 


( 8 ) 


(9) 


N <- 

Here, ipi is the relative phase of each transmitter, (y>(/3)) = 0 is the relative phase 
of the central transmitter (this is zero by construction as the central radiator 
is used as a reference) and /3 represents the angle corresponding to each of the 
points in which the focal plane was discretized (Fig. [2]). When setting Eq. © 
equal to zero, its analytical expression is a fourth-order periodic function of sin /3 
and cos/3 with coefficients depending on Xi , yt , and R\ (i = 1,... ,N e ). The 
apodization of the feed illumination limits the domain of interest to 0 < (3 < tt. 
In numerically exploring this function over this range, a maximum of four real 
roots could be identified, which repeat themselves with a periodicity of 27 t. 
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2.2. Optimization results 

The solution to the minimization problem defined in Eq. ([8]) was found with 
a quasi-Newton algorithm [M] and can be seen in Fig. [4] The emitters’ positions 
are indicated in red and present several characteristics similar to those found 
and discussed in [201 . First, they do not lie on the grating circle but on a curve 
that is tilted leftwards and intersects the grating circle at the central emitter 
before ending up inside the multimode region. Second, it was verified that 
the shape of this curve only approximates a circle with a radius ~ 2.2 R and 
is not symmetric. In the case presented here, this is caused by the absence 
of constraints on all stigmatic points. In general, for a two-stigmatic-point 
configuration the emitters lie exactly on the grating circle, as previously reported 
in the literature H13- The imposition of zero RMS phase error on a third 
stigmatic point (the blaze point) also caused a similar tilting effect [20], unless 
the emitters could be constrained to lie within a small distance (e.g., Ai/8) from 
the 2 R circle [19l . 



Figure 4: Optimized multimode region design with a resolving power TZ = 257 
and order M = 1. 
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The RMS phase error is shown in Fig. [5] Its values remain below 0.1 rad 
over an angular range spanning approximately from 16° to 176°. This phase 
error does not lead to a significant defocus of the light in the spectrometer [251 
and represents a ~ 30% improvement in the focal plane utilization over the 
previous designs [IS H3- Four stigmatic points are visible, but they are no 
longer associated with a predefined frequency, given the absence of constraints 
on them. In Fig. [5] it can be seen to what frequencies they correspond in this 
design. The nominal spectral range indicated in Table[2]is only partially covered 
down to ~ 510 GHz at 0°, while above 130° frequencies higher than 650 GHz 
show up, where the dispersion in the superconducting niobium’s reactance is no 
longer negligible |2B] . A slightly larger multimode region could be employed to 
reduce these tensions. 



Figure 5: RMS phase error distribution on the focal plane. The worst peak 
value is 0.1 rad at almost 45°, and four stigmatic points can be seen over a 160° 
angular range. 
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2.3. Power coupling efficiency 

The power coupling efficiency in the new design configuration was computed 
with the model described in [13 Section 4], The ratio of the power emitted by 
the feed horns to the power received by the antennas is approximately equal 
to unity. This high efficiency is the result of the absence of any higher-order 
diffraction peaks in the multimode region due to the relative emitter diameter, 
7 1 , being smaller than 1/2 [27] . 

The detailed coupling efficiency of the receiver feeds was not treated. For 
simplicity, the individual feed structures were mathematically modeled as if 
the focal surface were subdivided into apertures of equal size. The current 
configuration outperforms our previous first-order designs [13 ED] in terms of 
efficiency, while simultaneously providing the desired resolving power. 

3. Conclusions 

A design methodology was developed for high-resolution configurations of 
the /jl- Spec multimode region. The design procedure first maximizes the re¬ 
solving power subject to constraints on geometry, operating frequency range, 
and performance, thereby determining the order of the spectrometer. This then 
allows the RMS phase error on the instrument focal plane to be minimized. 
This work discussed a particular design achieved without constraining the RMS 
phase error to vanish at preselected points on the focal plane. This led to a 
configuration with four stigmatic points on the focal plane, a feature which can 
be used to increase the number of spectrometer channels as the phase error is 
reduced over a larger angular and spectral bandwidth. The design achieves a 
maximum RMS phase error equal to 0.1 rad, near-unity coupling, and a reso¬ 
lution of 257 in first order. Future work will be aimed at employing this design 
methodology to generate higher-resolution (JZ > 500) configurations. 
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